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Abstract 

In the first part of this review we introduce the basics theory behind 
geometric phases and emphasize their importance in quantum theory. The 
subject is presented in a general way so as to illustrate its wide applicabil- 
ity, but we also introduce a number of examples that will help the reader 
understand the basic issues involved. In the second part we show how 
to perform a universal quantum computation using only geometric effects 
appearing in quantum phases. It is then finally discussed how this geo- 
metric way of performing quantum gates can lead to a stable, large scale, 
intrinsically fault-tolerant quantum computer. 

1 Introduction 



This review is on the nature of geometric phases in quantum mechanics and 
their use in quantum computation. I will introduce the notion of a geometric 
phase in a general way which is neither traditional nor historical. This will, I 
believe, help us appreciate the importance of this notion and also its generality - 
it is a phenomenon that has been discovered and exploited to explain a number 
of effects in many different areas of physics (see the collection of papers in [1]). 
One of the potential applications of this phase which is of interest to us is in 
quantum computation (the standard text is [2]). By using geometric phases to 
implement any elementary quantum gate we can therefore perform any quan- 
tum computation. The main reason why we would want to use geometrical 
phases for quantum computation is that frequently geometrical evolutions are 
easier to control and may also be more resistant to random noise coming from 
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the environment. We will see, in Section 4, how any quantum computation can 
be implemented by using purely geometric (topological) phases. This review 
is nether complete nor detailed - I have tried to emphasize points and connec- 
tions that are usually neglected in literature, and to arrive at a potentially very 
important modern application in quantum computing. I have left a number of 
exercises and problems for the readers to address on their own. My intention is 
for the reader to acquire the basic flavor of the field which will then provide a 
good basis for a more independent study and work. 

The review is organized as follows: In section 2, I will review the notion 
of the geometric phase for the gauge invariant perspective. This section will 
also present a method of measuring geometric phases and discuss the classical 
analogue of the geometric phase. Various links with differential geometry will be 
emphasized. In section 3, I review some basic notions in quantum computation 
and present some simple algorithms. This is a necessary background to be 
able to show why geometrical phases can be used to implement any quantum 
computation. Section 4 then presents the most general way of implementing 
geometric quantum evolution, via the so called non-Abelian geometric phases. 
An example is given of a system that can support such an implementation. 
Section 5, finally, summarizes the review and offers an outlook on the field. 

2 Geometric Phases 

Quantum states are represented as vectors in a complex vector space (with a few 
other properties we need not worry about here). However, these vectors are only 
defined up to a unit modulus complex number, or a phase in other words. So, 
the two states |£) and e v |£) are indistinguishable as far as quantum mechanics 
is concerned - no set of measurements can discriminate them. It is interesting 
that although this statement looks fairly innocent at first sight, it can lead to 
some very profound conclusions. The route we will take here and that turns 
out to be very fruitful is to think of this extra phase in a geometric way as well 
- where we represent the phase as a unit length vector in the complex plane 
diagram. 

Although we think of the amplitudes between two quantum states as fun- 
damental entities in quantum theory, it is only the corresponding probabilities 
that we can ever observe experimentally Therefore, given two states \ipi) and 
\tpf), we can only measure the mod square of their overlap, \{4>i\ipj)\ 2 . What 
happens then if we wish to know the relative phase between these two states? 
Can we define their relative phase in spite of the fact that the absolute phase has 
no observable consequences? One idea to do so may be to look at the amplitude 
between the two states in the polar decomposition: 

(Vi|V/> = re ie *i (1) 
and define the phase as the relative phase between the two states. The main 
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Figure 1 : How do we compare phases of two different states when the probability 
between the states is the only experimentally well defined quantity in quantum 
mechanics? 



problem with this definition is that the states e ta \ipi) and e l/3 |-0/) which differ 
from the original states by an overall arbitrary phase will have a different relative 
phases by the amount of A8 — a — (3. This is not very satisfactory as there 
are infinitely many choices for A9 and they all look equally appropriate. We 
formally say that this definition of phase is gauge dependent (meaning "phase 
dependent") - a feature that is considered negative as all observable quantities 
should be independent of the choice of the gauge. How about defining a path 
connecting the two states, \ip(s)), such that when s = we have |V>i) and when 
s = 1 we have \ipf) (see Fig 1)? Then we can transport the states |V>( S )) from 
the position i to the position / and see how different the final phase is to that of 
\ipf) through interference. And if the states \ipf) and \ipi) transported to \ipf) 
interfere constructively, then they are in phase and if they interfere destructively, 
they arc out of phase; the degree of interference can therefore define the phase 
difference (I will be as vague as possible about what kind of interference this is 
- I will be much more precise later on in subsection 2,5). But which path do 
we choose and how do we know that the transport itself doesn't introduce any 
additional "twists and turns" in the phases so that we are actually comparing 
some different phases to the original ones? This, in fact, is a well known problem 
in (differential) geometry (sec [3] for an excellent introduction). Let's state the 
problem in an abstract setting first (those of you familiar with general relativity 
will know immediately what I am taking about) and then specialize to quantum 
geometric phases. 

2.1 Parallel Transport 

Suppose we have a curved manifold (as in Fig. 2) and we have a vector at a 
point A and another at a point B. What is the angle between the two vectors 
(since phases in quantum mechanics can be thought of as vectors, this question 
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Figure 2: Comparing direction of two vectors at two different points is a non- 
trivial problem on a curved manifold. We have to somehow transport the vectors 
to the same point and then see what the angle between them is. Otherwise, 
comparing two vectors at different locations is difficult. 



is the same as our original question of relative phase between two quantum 
states). So, how do we measure the angle when the two vectors are at different 
positions? Well, we can transport one of them to the other one and then, when 
they are next to each other, the angle is easily measured (in the usual way 
as the "angular distance" between the vectors). But, again, we don't want 
the transport itself to introduce any additional angles - we'd like it to be as 
"straight" as possible. The straightest possible path is known as a geodesic and 
the corresponding evolution along this path is known as the parallel transport 
(the usefulness of this concept in physics cannot be overemphasized - see [4] ) . To 
define a parallel transport, let's look at the infinitesimal evolution, from \tp(s)) 
to \ip(s + ds)}. If we don't want there to be any twists and turns in the phase, 
even infinitesimally, then the two states should be in phase. So, we require that 

Arg{<V>(s)|V(s + ds)>} = (2) 

This is the same as asking that (tp(s)\ip(s + ds)} is purely real, i.e. 

Im{^(s)|V(s + da))} = Im{<V>( S )MV(s)}} = (3) 

up to second order. But {ip(s)\d\'ilj(s)) is purely imaginary anyway (prove it by 
differentiating (ip(s)\ip(s)) = 1 by ds), and hence the parallel transport condition 
becomes: 

(V(s)W(s)) = (4) 
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If the evolution satisfies this equation, then the phase is said to be parallcly 
transported, which is what we mentioned was required in order to be able to 
define a relative phase between two states (also see [5] for an accessible treatment 
of the phase parallel transport) . This definition of parallel transport is, however, 
not automatically gauge invariant. By this we mean that if instead of the 
state \ip(s)), we use the state |f/>(s)) = e la ^\tp(s)) , then the parallel transport 
condition changes by the amount 

${s)\d\${8)) = W(s)\d\1>{8))+i^ds (5) 

as can easily be checked. In order to obtain something that is gauge invariant we 
can integrate the expression (ip\d\ip) over a closed loop, giving us the expression 
for the geometric phase, and then exponentiate the result (the integral of ^ 
over a closed loop is 2ir, whose exponential is equal to unity). So, the geometric 
phase resulting from the parallel transport is 

7 = j\^{s)\j- s \^{s))ds (6) 

and its exponential (over a closed loop) is gauge independent, but not path in- 
dependent . Incidentally, note that we are requiring that the phase difference 
between every two infinitesimal points is zero. So, how come, if every infinites- 
imal difference is zero, that on a finite stretch a non-zero phase builds up? The 
answer is that the underlying space is curved and it is the curvature that is 
reflected in the phase difference; in fact, the curvature is equal the phase differ- 
ence (up to a constant factor) . When a quantity vanishes infinitesimally, but its 
integral over a finite region doesn't, then this quantity is called non-integrable. 
Therefore, it can be succinctly said that geometric phases are a manifestation 
of non-integrable phase factors in quantum mechanics. Let's look at a spin half 
system to illustrate this point. 



2.2 The Bloch Sphere 

Two level systems are ubiquitous in nature and can be used to store and ma- 
nipulate any form of information (in quantum computing they are known as 
quantum bits or qubits). Two level systems are very conveniently represented 
on a sphere, such that all pure states lie on the surface of this sphere and all 
the mixed states are inside. There is a one-to-one correspondence between the 
points on and inside the sphere and quantum states of a two level system (Fig 
3). This is because we need at most three parameters to specify every two level 
state, as can be seen from the general density matrix for the systems: 

i 

x Path dependence cannot easily be eliminated, although this is the key idea behind con- 
structing "topological" as opposed to "geometric" phases 
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where a are the Pauli spin matrices and Si = Xxoip are the x, y, z coordinates 
on the sphere. If a state is pure than s\ + s 2 y + si = 1 which, as we said, is a 
point on the surface of the sphere (a pure state in other words) . 



Figure 3: Parallel transport of the phase on the Bloch sphere. The presented 
evolution can be implemented in many different ways, one of which is the adia- 
batic Schrodinger equation as explained in the next section. The phase vector, 
initially at the North Pole and lying in the plane of the paper, for example, 
is then transported to the equator. After that, it moves on the equator, and 
makes a 90 degrees shift, before coming back to the North Pole. The final phase 
vector, sticking out of the page, is at an angle to the initial vector which is also 
90 degrees. 

Suppose that we now evolve from the state |0) to the state |+) = |0) + |1), 
then to | a;) = |0) + and finally back to the state |0). On the Bloch sphere 
we are going from the North Pole to the Equator, then we move on the Equator 
by an angle of ir/2 and finally we move back to the North Pole). What is the 
corresponding geometric phase? To see this we start with a tangential vector 
initially at the North Pole pointing in some direction (there are infinitely many 
directions corresponding to infinitely many arbitrary phases to start with). If 
we now parallely transport this vector along the described path, then we end up 
with a vector pointing in a different direction to the original one (even though, 
infinitcsimally, the phase vector has always stayed parallel to itself!). The angle 
between the initial and the final vector is n/2, which is exactly equal to the area 




Bloch Sphere 
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covered by the state vector during the transport, or the corresponding solid 
angle of the transport (the two are just different ways of phrasing the same 
thing and can be seen to be equivalent using Stokes' theorem). Therefore, the 
two orthogonal states |0) and |1) evolve in the following way: 

|0) -» e^/ 2 |0) (8) 
|1) -» e- jf2 / 2 |l> (9) 

where f2 is the solid angle (the factor of half is there because the orthogonal 
states are n away from each other in the Bloch representation). This is im- 
portant as it shows that orthogonal states pick up opposite phases of equal 
magnitude. The phase can be computed from the formula in eq. (6), but we 
note that it can also be written in a "discrete" way as 

arg{(0|+)(+| a; )( a ;|0)} (10) 

This beautiful formula - originally due to Pancharatnam [1] - will be the basis 
of our general formulation of the geometric phase for pure states. From it we 
will be able to find that the space of quantum two level systems is curved. Are 
there, for that matter, any quantum states whose space is not curved (i.e. is 
flat) and how do we determine the curvature in general? If you want to measure 
a curvature at a point, then make a small loop around that point and compare 
the initial phase of your state to the phase after completing this loop. Suppose 
that the evolution is \ip(si)) — > \ip(si + 5 s)) — ► \tp(si + 5s + As)) — > \ip(si)). 
Then, the final phase difference is (according to Pancharatnam) 

50 = arg{(V(si)|V(si + Ss))(^( Si + 5s)\^( Si +Ss + Ss'))(^( Si + Ss + Ss')\^( Si ))} 
= arg{l + i(A Wa .)|^.|^.))_^. Ws .)|^|^( a .)))} (11) 

= ^M'OI^IVM-^WOI^M*)) (12) 

The definition of curvature, K, is 50 = K5A, where A is the area enclosed (this 
definition of curvature is due to Gauss). You can convince yourselves that for 
spin half particles SO = 5 A, which means that the curvature is 1 (it is really 1/r 2 
as for any sphere, but the Bloch sphere has a unit radius). For optical coherent 
states, on the other hand, SO = (if we vary the phase and the amplitude of the 
state only), indicating that they line on a flat surface (which you can prove as 
an exercise, but see [6] for more details on the Riemannian structure of quantum 
states). 

This formalism is very powerful and important in the modern formulation 
of physical theories (the so called gauge theories - see [4] for more details). But 
before we admire the beauty and cleverness of all this any further, let's ask 
ourselves a very basic question: how do we actually (physically) implement the 
parallel transport given above? 
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2.3 Adiabatic implementation of parallel transport 

Here we present one such implementation due to Berry [7], who actually dis- 
covered the quantum geometric phase in 1984 in this way. A wavefunction of 
a system is a function of a set of parameters (which generally depend on time) 
and the time itself, |\&(s(i), t)). Suppose that the Hamiltonian is a function of 
s only H = H(s{t)). Suppose, in addition, that these parameters change very 
slowly so that the system, which starts in a eigenstate of the initial Hamiltonian, 
stays an eigenstate of the instantaneous Hamiltonian, i. e. 

H(s(t))\*»(s(t), t)) = E n (s(t))\* n (s(t), t)) 

(this is what is known as the adiabatic approximation in quantum physics- 
the other extreme is the sudden approximation where the Hamiltonian changes 
instantaneously). The Schrodingcr equation for this system is 

if t \* n (s(t),t))=H(s(t))\* n (s(t),t)) 

(where we assume throughout that K = 1). We would like to show that the 
adiabatic evolution naturally implements the parallel transport of the phase 
of the quantum state. Note that apart from the adiabaticity assumption our 
analysis is completely general. Multiplying the Schrodinger equation by (^ n \, 
and taking the eigenvalue equation into account we obtain: 

*<*n|^|*n) =E n 

Now, every state gains a dynamical phase (as well as a geometric one) as it 
evolves. We'd like to get rid of the dynamical component of the phase and not 
take it into account when we are measuring the geometric contribution. To do 
so we define a new wavefunction by taking into account the dynamical phase 

|$„(«(t), t)> := exp{iE n (s(t))}\* n (s(t), t)) 

which satisfies the following equation: 

<$(s(i),t)|||*(s(t),t)>=0 

And this is the same as the parallel transport condition. So, the geometric part 
of the quantum phase is parallely transported when the state evolves under the 
Schrodingcr equation in the adiabatic approximation. So can we derive a closed 
formula for the geometric phase? The answer is yes, and the solution to this 
equation is given by: 

where 



8 



We notice that the right hand side involves the 1-form (3 describing the (symplec- 
tic) manifold of the projective Hilbert space 2 - see Simon's original differential 
geometric interpretation [8]. Therefore, 



where (3 = — i($(s(t))\j^\$(s(t))). Integrating over a closed circuit OS we solve 
for the Berry phase 

7 = / P 

JdS 

Invoking Stokes' theorem this becomes 




where S is the area enclosed by the circuit dS and a is the 2-form on the sym- 
plcctic manifold of the projective Hilbert space. This phase is called the Berry 
phase (the Berry and geometric phase are terms that I will use interchangeably, 
although this is by no means the standard convention in the literature). 

In summary, the quantum phase arising in a general evolution of a quantum 
system consists of two parts: the dynamical part given by the formula 




and the geometric part, given by 

7 = /<*(-)!>(<)>!* 

The total phase is then the sum of the two contributions, <fi = S + 7. 

We can encapsulate this whole subsection by saying that 

The Schrodinger equation in the adiabatic limit implements the parallel trans-% 
port of the phase of the evolving state. 

In order to become more acquainted with the notion of geometric phases, 
we look at the classical analogue of this phenomenon. 



2.4 Classical geometric phase 

In this subsection we will talk about astronauts, cats and the Foucault pendulum 
in this section, all of which are classical objects. What do all these have in 

differential forms are just things that exist under the integral sign, however, in mathe- 
matics they can be given a life of their own. We will not need this formalization in the rest 
of this review, but they feature very strongly in modern theoretical physics, which is why I 
mention them here. 
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common? They all have the geometric phase at the root of their basic behaviour! 
Suppose an astronaut is in free space and initially facing away from his space 
ship. He wants to turn around by 180 degrees to return to the ship, but there 
is nothing around him to push against or hold onto in order to turn around. 
His task may at first sight look impossible to achieve: his angular momentum 
is initially zero and in order to turn he seems to need to generate some angular 
momentum. But there arc no forces externally, so this looks indeed impossible! 
However, anyone who has been observing cats knows that this reasoning is false. 
Cats face the same problem everyday and they manage to solve it. If a cat falls 
from a certain height and was initially upside down, she is still going to (almost 
always) land on her feet (hence the allegorical expression: "cats always land on 
their feet"). So, how do cats manage to apparently "violate the conservation of 
angular momentum" ? Well, they don't violate angular momentum conservation: 
they don't use dynamics to turn, they use topology! Fig. 4 shows how the 
astronaut does something similar to the cat to turn around. Note that his hand 
movement is the same as the spin half evolution on the Bloch sphere - hence 
the angle by which he turns can be calculated in the same way. 



Figure 4: A person rotating while standing on the chair by using a purely 
geometrical effect while maintaining a zero angular momentum. In the position 
three, the person turns his arms around, and at the same time his body has 
to turn in the opposite direction to maintain zero angular momentum on the 
whole. Therefore by repeating the whole cycle a number of times, the person 
can make a rotation on the chair by any desired angle (try it for yourself!). This 
same effect is used by a cat to turn around in the air and land on the feel, even 
when initially being upside down. 

Interestingly enough, the Foucault pendulum can be explained in the same 
geometric way [9]. Imagine that a pendulum is suspended at some latitude 9, 
and that we observe its swinging motion as the Earth spins around its own axis. 
It is well known that the pendulum will after one rotation of the Earth swing in 
a plane at an angle to the original plane of motion. To see this, let us obtain an 
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equation of motion of the pendulum. The Lagrangian for this problem is given 
by 

1 dx 2 dy 2 1 2 2 2 
L= 2 m( ^ + ^ ) -2 m " (a; +2/ )-^cos^-- y -) (13) 

where f2 is the Earth's frequency of rotation (2tt per day) , m is the mass of the 
pendulum and u its natural frequency of swinging. The corresponding Euler- 
Lagrange equations of motion are easy to obtain (exercise!) and their solution 
written in the coordinate z = x + iy is (in the adiabatic limit): 

z(t) wi e" inOTW e- iwt (14) 

This form of the solution is convenient as it allows us to see two phases in the 
motion of the pendulum: the dynamical phase ut which is the same as the 
quantum dynamical phase seen before and the geometric phase f2 cos 9t. Note 
that after one day revolution of the earth, the geometric phase is 2n(l — cos 9) 
which is the same as the spin half geometrical phase. This is because our planet 
is a very good, albeit somewhat imperfect, physical representation of the Bloch 
sphere! 



2.5 Experimental realization: The Pancharatnam connec- 
tion 

I now summarize the theory so far: we would like to compare phases of a quan- 
tum state at two different points. To do so we have to chose a path connecting 
the two points, and evolve one of the points along the path until it coincides 
with the other one. Then we interfere them to infer the phase. But how does 
this interference work? 

Consider a conventional Mach-Zehnder interferometer in which the beam- 
pair spans a two dimensional Hilbert space Tri = {|0), |1)} [10]. The state vectors 
|0) and |1) can be taken as (orthogonal) wave packets that move in two given 
directions defined by the geometry of the interferometer. In this basis, we may 
represent mirrors, beam-splitters and relative U(l) phase shifts by the unitary 
operators 

Um = (J J ) , U B = i= ( \ _\ 
U(l) = ( ^ J ) , (15) 

respectively. An input pure state pi n = |0)(0| of the interferometer transforms 
into the output state 

Pout = U B U M U{l)UBP l nUlu\l)Ul I ul 

1 / 1 + cos x i sin x 

2 \ — isinx 1 — cosx 



(16) 
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that yields the intensity along |0) as I oc 1 + cosx- Thus the relative U{1) phase 
X could be observed in the output signal of the interferometer. 

Now assume that the particles carry additional internal degrees of freedom, 
e.g., spin (it is this degree of freedom that will be parallely transported). This 
internal spin space Hi = C N is spanned by the vectors |fc), k = 1,2,. ..N. The 
density operator could be made to change inside the interferometer 



Po 



U iPo U} 



(17) 



with Ui a unitary transformation acting only on the internal degrees of freedom 
(we will see later on that this transformation need not be unitary, but could 
be a more general completely positive map). Mirrors and beam-splitters are 
assumed to leave the internal state unchanged so that we may replace Um and 
U B by U M = Um ® 1« an d U B = Ub <£> T, respectively, lj being the internal 
unit operator. Furthermore, we introduce the unitary transformation 



U = 




1 



■u l + 







h, 



(18) 



The operators Um, Ub, and U act on the full Hilbert space H <g> Hi- U 
corresponds to the application of Ui along the |1) path and the U(l) phase \ 
similarly along |0). We shall use U to generalize the notion of phase to unitarily 
evolving mixed states. 

Let an incoming state given by the density matrix g in — p in ® p — \ 0) (0 1 ® po 
be split coherently by a beam-splitter and re-combine at a second beam-splitter 
after being reflected by two mirrors. Suppose that U is applied between the 
first beam-splitter and the mirror pair. The incoming state transforms into the 
output state 

g out = U B U M UU B ft n Ut,Utukl4. (19) 



Inserting Eqs. (15) and (18) into Eq. (19) yields 



Qout 



+e 



1 -1 
1 -1 



U lPo Uj + 
®PoU} 
<g> Uipo 



1 

-1 



'Po 



(20) 



The output intensity along |0) is 

I oc Tr(Uip Uj + Po + e- ix U i po + e ix PoUfy 
oc l + \Tr {Uipo) | cos [x - arg Tr (U lPo )} , 

where we have used Tr(p uj) = [Tr (Uip )]* ■ 



(21) 
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X 



7 ^-d^ 1 



|0). 



7^ 



Ui 



Figure 5: The Mach Zehnder interferometer intended to implement the mea- 
surement of the Pancharatnam phase. The internal degree of the system going 
through the interferometer is manipulated only in one branch of the interfer- 
ometer. The manipulated branch is then interfered at the end with the non- 
manipulated branch to infer the phase difference. 



The important observation from Eq. (21) is that the interference oscillations 
produced by the variable U(l) phase \ * s shifted by 4> = argTr (UiPo) for any 
internal input state p , be it mixed or pure. This phase difference reduces for 
pure states po = \ipo)(ipo\ to the Pancharatnam phase difference between Ui\tpo) 
and \ipo). Moreover the visibility of the interference pattern is v = \Tr {UiPo) \ > 
0, which reduces to the expected v = |(V>o|£^|V'o)l for pure states. 

It is clear from this experiment what it means to have a parallel transport 
of phase, i.e. no phase introduced during an infinitesimal evolution. We require 
that 

argTr (U(t, t + dt)p Q ) = (22) 

which is the same as the parallel transport condition presented before in cq. 
(4). Suppose for that matter that we have no dynamics and only projective 
transformations (the transformation U is in this case no longer unitary, but the 
same above analysis still holds). Then U — |^> n )(^ n _i|...|^2)(V'i| so tnat the 
formula for the phase now becomes 

7 = ar<7{<V>i|V>2}<V>2|V>3> ■ • ■ (V>n|V-i)} (23) 

which is known as the Pancharatnam phase. Now, this is the most general 
form of the phase for pure state, which leaves the motion between the states 
completely arbitrary: it could be adiabatic, and cyclic, but it could also be 
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non-unitary and open. There are many advantages of expressing the phase in 
this way: 

• it is aesthetically pleasing! 

• it is manifestly gauge invariant since whenever we have the ket \ipj) in the 
expansion, we also have its bra present, meaning that all the extra 
phases of individual states cancel. Thus, states only appear as projectors 
which is the key to achieving independence of phase. As a slight digression, 
those of you who know a bit about the lattice gauge theory (see [11] for 
an excellent succinct introduction) will immediately notice the similarity 
with Wilson's loop construction. In order to obtain the action leading to 
the geometric phase we have to sum up the infinitesimal loops, Wp, over 
all points P. Each loop is defined as 

W P = (i/>{xp)\i/)(xp + dx))...{ip(x P +dx + 5x)\tp{x P )) (24) 

After a simple manipulation and using the Baker-Hausdorff formula 

exp(AA) exp(AS) = exp(A(A + B) + A 2 /2[A, B] + 0(A 3 ), 

it can be shown that 

W P = e iF ^ dx ^ 5x " (25) 
where the fictitious geometric field strength is 

F/J.V = df J/ A v — d v A^ — i[A mu , A nu ] (26) 

and the geometric potential (also called connection in the lectures) is 

A v = ty{x)\d v \1>{x)) (27) 

To obtain the Wilson loop over a finite range we just need to add up all 
the infinitesimal loops such that W — ~^2 P Wp. Pushing the analogy with 
gauge theories we can think of a fictitious topological charge Q T defined 
via Gauss' law: 

j j> B^dA = Q T (28) 
Also, a covariant derivative of a field $ is given by 

= d a <$> + iA a <$> (29) 
and the field satisfies the following well known identity 

Fp.v\p + Fvp;n + Fvp\n = (30) 
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• the formulation is very general. Any most general quantum evolution 
(the so called Completely positive, trace preserving map) going through 
the states above will lead to the same phase. There are many issues in this 
generalization that are still unclear, but I will not make any comments on 
this as it would obscure the exposition (see Sjoqvist et al. [10] for details). 

• offers us a way of generalizing the phase to mixed states (as averaging over 
the corresponding pure states). 

So, the geometric phase is independent of dynamics and gauge-free, however, 
it does depend on the underlying geometry of the evolution. To what extent 
can we eliminate the sensitivity to geometry? 

2.6 Topology: The Aharonov-Bohm effect 

The Aharonov-Bohm effect [12] can be revisited and explained in terms of the 
geometrical phase notion. Whereas the curvature in the above examples was in 
some sense fictitious (i.e. not generated by a "real" field, but by the intrinsic 
curvature in the geometry of quantum states), the Aharanov-Bohm effect is a 
manifestation of the electro-magnetic field. From this perspective we can even 
say that 

electro-magnetism is the gauge invariant manifestation of a non-integrable 
quantum phase factor. 

In other words, this means that the electromegnetic field (more precisely 
the potential) can be derived by postulating that the Schrodingcr equation is 
invariant under local phase change of the wavefunction. 

The Aharonov-Bohm set-up is the best example that explores the effect of 
the electro-magnetic field on phase through a beam of electrons in a double 
connected region where the field strength is zero as in Fig. 6. The phase gained 
by the beam is 

r — ie$, .„„. 
exp{ — } (31) 

where 3> is the flux of the field. This phase has been experimentally observed. 
We see that the phase is path dependent (which is another way of saying that 
it is non-integrable). Even though it is path dependent, if we make one round 
trip about the flux, then we always gain the same amount of phase no matter 
what the path is. This a nice topological feature of the Aharonov-Bohm phase; 
the phase is then dependent on the number of times the flux has been encircled 
(the so called winding number in topology). It is this topological feature that 
we would like to have when we implement quantum computation (this idea was 
first presented in a seminal paper by Kitaev [13]). However, we would like the 
topology to be of a non-abelian (i.e. non-commutative) character as we explain 
later on in the review. 
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Figure 6: The Aharonov Bohm set up. A beam of electrons undergoes a double 
slit experiment with a solenoid placed behind the slits. The soleniod generates 
a magnetic field which is confined to the small region of the solenoid and is 
responsible for a shift in the interference pattern at a screen placed some large 
distance away. 



3 Quantum Computation 

So far we have explored connections between physics and geometry through the 
concept of the geometric phase which we saw is manifested both in quantum 
and classical mechanics. Now we plunge into another equally exciting and cross- 
disciplinary area: quantum information processing (see [2] for details). 

All information is ultimately encoded in the most elementary quantum sys- 
tems, so called quantum bits. A qubit is a quantum system in which the Boolean 
states and 1 are represented by a prescribed pair of normalised and mutually 
orthogonal quantum states labeled as {|0), |1)}. Unlike a simple Boolean vari- 
able, a qubit, typically a microscopic system such as an atom, a nuclear spin, or 
a polarised photon, can exist in an arbitrary superposition a\0) + (3\1), making 
it more powerful as a computational resource. 

In quantum computation, we set some register of qubits to an "input" state, 
evolve the qubits unitarily using simple building-block operations and then take 
the final state as "output" . More formally, a quantum logic gate is a device 
which performs a fixed unitary operation on selected qubits in a fixed period 
of time and a quantum network is a device consisting of quantum logic gates 
whose computational steps are synchronised in time. The outputs of some of 
the gates are connected by wires to the inputs of others. Also we define the size 
of the network to be the number of gates it contains. 

3.1 Quantum logic gates and universality 

The most common quantum gate is the Hadamard gate, a single qubit gate H 
performing the unitary transformation known as the Hadamard transform. It 



1G 



is defined as 



H = 



1 

V2 



x) — H — 



(-if\\x)+\i-x) (32) 



The matrix is written in the computational basis {|0), |1)} and the diagram on 
the right provides a schematic representation of the gate H acting on a qubit 
in state \x), with x = 0, 1. 

The addition of another single qubit gate, the phase shift gate (j), defined as 
| 0) h-> | 0) and 1 1) e 1 ^ | 1), or, in matrix notation, 



1 





\x) 



a ix<p 



X) 



(33) 



is actually sufficient to construct the following network (of size four), which 
generates the most general pure state of a single qubit (up to a global phase), 



10} 



H 



29 



H 



cos0|O) +e i0 sin0|l). 



(34) 



Consequently, the Hadamard and phase gates are sufficient to construct any 
unitary operation on a single qubit. We have used both the above gates in our 
discussion of the Pancharatnam phase in the Mach-Zehnder interferometer. 

Thus the Hadamard gates and the phase gates can be used to transform the 
input state |0)|0)...|0) of n qubits into any state of the type \^2}--- \^n), 
where is an arbitrary superposition of |0) and |1). These are rather special 
n-qubit states, called the product states or the separable states. In general, a 
register of n qubits can be prepared in states which are not separable, known 
as entangled states. 

However, in order to entangle two (or more qubits) it is necessary to have 
access to two-qubit gates. One such gate is the controlled phase shift gate B{<j>) 
defined as 



B{4>) 

















1 














1 





V o 











\x) 

\y) 



} e lxv4, \x)\y). (35) 



The matrix is written in the computational basis {|00), |01), |10), |11)} and the 
diagram on the right shows the structure of the gate. Note: we know how to 
implement the single qubit phase gate geometrically, but not the other two gates 
aforementioned. Which gates are important to be able to implement in order 
to perform any quantum computation? 

An important result in the theory of quantum computation states that the 
Hadamard gate, and all B(<j)) controlled phase gates form an universal set of 
gates: if the Hadamard gate as well as all B((f>) gates are available then any 
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n-qubit unitary operation can be simulated exactly with less than CA n n such 
gates, for some constant C [2]. Consequently, being able to implement 1- and 
2- qubit phase gates is of crucial importance in quantum computation. In the 
remaining part of the lectures we describe a new method for implementing the 
controlled phase gates based explicitly on geometric phases rather than dynamic 
ones. Thus if we show that the universal set can be implemented geometrically, 
then any computation can be implemented geometrically (see [14] for an NMR 
geometric implementation of the controlled phase gate). 

3.2 Lightning review of quantum algorithms 

Deutsch's problem [15] is the simplest possible example which illustrates the ad- 
vantages of quantum computation [16]. The problem is the following. Suppose 
that we are given a binary function of a binary variable / : {0, 1} — ► {0, 1}. 
Thus, /(0) can either be or 1, and /(l) likewise can either be or 1, giving 
altogether four possibilities. However, suppose that we are not interested in the 
particular values of the function at and 1, but we need to know whether the 
function is: 1) constant, i.e. /(0) = /(l), or 2) varying, i.e. /(0) ^ /(l)- Now 
Deutsch poses the following task: by computing / only once determine whether 
it is constant or varying. This kind of problem is generally referred to as a 
promise algorithm, because one property out of a certain number of properties 
is initially promised to hold, and our task is to determine computationally which 
one holds. 

First of all, classically finding out in one step whether a function is constant 
or varying is clearly impossible. We need to compute /(0) and then compute 
/(l) in order to compare them. There is no way out of this double evalua- 
tion. Quantum mechanically, however, there is a simple method to achieve this 
task by computing / only once! We need to implement the following unitary 
transformation 

\x) => e l7rf(x) \x) 

where x — 0,1. Now we can see where the power of quantum computers is 
fully realised: each of the components in the superposition of |0) + |1) undergoes 
the same above evolution "simultaneously", leading to the powerful "quantum 
parallelism". This feature is true for quantum computation in general. Let us 
look at the two possibilities now: 

1. if / is constant, then the final state is 

|*out>=±(|0> + |l». 

and 

2. if, on the other hand, / is varying, then, the final state is 

|*out>=±(|0>-|l». 
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Note that the output qubit emerges in two different orthogonal states, depending 
on the type of /. These two states can be distinguished with 100 percent 
efficiency. This is easy to see if we first perform a Hadamard transformation 
on this qubit, leading to the state |0) if the function is constant, and to the 
state |1) if the function is varying. Now a single projective measurement in 
0, 1 basis determines the type of the function. Therefore, unlike their classical 
counterparts, quantum computers can solve Deutsch's problem. 

It should be emphasised that this quantum computation, although extremely 
simple, contains all the main features of successful quantum algorithms: it can 
be shown that all quantum computations are just more complicated variations of 
Deutsch's problem. In particular, Shor's algorithm for factorisation of numbers 
is just a scaled up version of Dcutsch as it is based on period finding. Note, 
finally, that the main operation in Deutsch's algorithm was the introduction of 
a phase conditional on the state of the qubit. This phase can, of course, be 
introduced in a purely geometrical way, by a spin-half encircling an area of the 
size of 2n. 

4 Topological implementation of quantum gates 

Speaking somewhat loosely, we would now like to construct a universal set 
of quantum gates, which we know is impossible if we only have the (abelian) 
geometric phases (because we cannot implement a Hadamard transform with 
them). However, in addition to achieving the universality, we'd like to be able 
to achieve a certain degree of resistance to errors, but without performing any 
error correction. How far can this idea be pushed? 

4.1 Non-abelian formalism 

The best way of explaining the non-abelian phases is to look at how they arise 
during quantum evolution. The term "non-abelian" refers to the fact that two 
different phases don't commute. Note that the ordinary (Berry) phases we 
discussed so far are abelian, i.e. 



This inequality, however, would not be satisfied if we had non-abelian phases. 
This therefore implies that the non-abelian phases would not be represented 
by ordinary numbers, but by matrices. We now show that this scenario arises 
during the Schrodinger evolution in the adiabatic setting when our system is 
degenerate. Suppose that we have 



e' 



H(t)\Mt)) 
H(t)\Mt)) 



E(t)\Mt)) 
E(t)\Mt)) 



(36) 
(37) 
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and (i^i(t)\tl>2(t)) = 0. From the Schrodingcr equation we can conclude: 

= Efaty) (38) 

i(H^\1>) = E(W1>) (39) 

We recast this pair of equations by writing \ip) — ai(t)\ipi) + a2(t)\ip2) , and 
subtracting the dynamical phases by putting E(t) = (for all time t) we obtain 

^+ai(Vi|^|Vi)+a2(Vi|^l^2) - (40) 
^+a 2 (Hj t \^)+a 1 (^\f t \^ 1 ) = (41) 
We can rewrite this equation as a vector-matrix equation, 



where 



v(t) = 



ai(t) 
a 2 (t) 



The formal solution to this equation can be written in the usual "Dyson per- 
turbative expansion" : 



v(s) ="Pexp{i / A 1 (s^dsjviO) 

J i 



where V indicated the path-ordering operation of non-commuting factors (s)Tj J 
This non-abelian phases is in general very difficult to evaluate, because of the 
path ordering operation. Note that if the state starts in the degenerate subspace 
it always remains there according to the adiabatic theorem 3 . 

I would like to mention potential advantages of using topological evolution 
to implement quantum gates. First of all, there is no dynamical phase in the 
evolution. This is because we are using degenerate states to encode information 
so that the dynamical phase is the same for both states (and it factors out as it 
were). Also, all the errors stemming from the dynamical phase are automatically 
eliminated. Secondly, the states being degenerate do not suffer from any bit flip 
errors between the states (like the spontaneous emission). So, the evolution is 
protected against these errors as well. Thirdly, the size of the error depends on 

3 We mention that the SU(2) Yang - Mills formalism is based on the non-abelian gauge 
invariance and can reproduce the weak interaction and the SU (3) formalism due to Gell-Mann 
can be used to represent the strong interaction. These are just generalizations of the U(l) 
formalism for Elactro-magnetism based of Weyl's ideas. The interested reader is advised to 
read [4] 
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the area covered and is therefore immune to random noise (at least in the first 
order) in the driving of the evolution. This is because the area is preserved under 
such a noise as formally proven by Palma [17]. Also, by tuning the parameters 
of the driving field it may be possible to make the phase independent of the area 
to a large extent and make it dependent only on a singular topological feature - 
such as in the Aharonov-Bohm effect where the flux can be confined to a small 
area - and this would then make the phase resistant under very general errors. 

So, in order to see how this works in practice we take an atomic system 
as our model implementing the non-abelian evolution. We'll see that quantum 
computation can easily be implemented in this way. The question, of course, is 
the one about the ultimate benefits of this implementation. Although there are 
some obvious benefits, as listed above, there are also some serious shortcomings, 
and so the jury is still out on this issue. 



4.2 Example 



Let's look at the following 4 level system analyzed by Unanyan, Shore and 
Bergmann [18]. They considered a four level system with three degenerate 
levels 1,3,4 and one level 2 with a different energy as in Fig 7. This system 
stores one bit of information in the levels 1 and 2 (hence there is double the 
redundancy in the encoding of information) . We have the following Hamiltonian 



Hit) 



where P,Q,S are some functions of time (to be chosen at leisure). It is not 
difficult to find eigenvalues and eigenvectors of this matrix (exercise!). There 
are two degenerate eigenvectors (with the corresponding zero eigenvalue for all 
times) which will be implementing our qubit and they are 
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P(t) 








\ 




p(t) 





S(t) 


Qit) 







S(t) 










V 





Q(t) 











V 



cos0(i) 


-sin 0(f) 




/ 



and 



$ 2 (t) 



/ sin (j){t) sin 6{t) 


sin 4>(t) cos 6{t) 
\ — cos 4>(t) 



where tan0(t) = P(t)/Q(t) and tan0(t) = Q(t)/^P(t) 2 + Q(t) 2 . 

In the adiabatic limit (which we will assume here) these are the only impor- 
tant ones as the state will remain in their subspace throughout the evolution 
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Figure 7: The four level system that can be used for non-abelian quantum 
computation to encode one qubit of information in two degenerate levels. The 
method is detailed in the text. 



(this is guaranteed by the adiabatic theorem). Although, in general, the Dyson 
equation is difficult to solve, in this special example we can write down a closed 
form expression [18]. The unitary matrix representing the geometrical evolution 
of the degenerate states is 



B{t) 

where 



cos 7(f) sin 7(f) 
— sin -f(t) cos 7(f) 



7 (f) = J sin0(f)^df (42) 

This therefore allows us to calculate the non-abelian phase for any closed path 
in the parametric space. After some time we suppose that the parameters return 
to their original value. So, at the end of the interaction we have the matrix 



B(oo) 

where 



cos 7/ sin 7/ 
— sin 7/ cos 7/ 



7 f = / Q (SdP - PdS) (43) 

s To {p 2 + s 2 yQ 2 + p 2 + s* K 

which can be evaluated using Stokes' theorem (the phase will in general de- 
pend on the path, as explained before). So, we can have a non-abelian phase 
implementing a Hadamard gate. With two systems of this type (mutually in- 
teracting) we can implement a controlled Not gate and therefore (at least in 
principle) have a universal quantum computer (see [19]). As an exercise I leave 
it to the reader to implement Deutsch's algorithm in the purely topological way. 



5 Conclusions and Outlook 

We have seen what geometric phases are and appreciated their importance in 
physics in general. They were shown to be linked to the modern ideas of gauge 
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theories and differential geometry, in particular the idea of parallel transport. 
In addition to their fundamental value, they also have a potential to implement 
quantum gates that are, at least in principle, more reliable than any existing 
gates. In order to do so I have argued that a high degree of topological de- 
pendence has to be achieved - i.e. independence of global geometry, which is 
certainly not that easy in practice. These topological gates will most likely have 
to be combined with other error correcting mechanisms in order to make a reli- 
able quantum computer. The exact engineering design will, in the end, depend 
on the actual medium that is chosen to support quantum computation. It is fair 
to say that when this review was completed (in December 2002) there was no 
fully fault tolerant implementation of quantum computation. Fortunately there 
has been a huge number of proposals to implement such computations with 
various solid-state systems (such as Josephson Junctions and Quantum Dots), 
Nuclear Magnetic Resonance, Trapped Atoms and so on, and the whole field is, 
therefore, very much alive at present. 
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